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Abstract 

We address the recent anomalies in semi-leptonic _B-meson decays using a model of fermion 
masses based on the 17(2) flavor symmetry. The new contributions to 6 —>■ si£ transitions 
arise due to a tree-level exchange of a Z' vector boson gauging a [/(I) subgroup of the flavor 
symmetry. They are controlled by a single parameter and are approximately aligned to the 
Standard Model prediction, with constructive interference in the e-channel and destructive 
interference in the /i-channel. The current experimental data on semi-leptonic i?-meson 
decays can be very well reproduced without violating existing constraints from flavor 
violation in the quark and lepton sectors. Our model will be tested by new measurements 
of & —>■ s££ transitions and also by future electroweak precision tests, direct Z' searches, 
and fi-e conversion in nuclei. 



1 Introduction 


While direct searches for new particles at the LHC have so far been inconclusive, recent 
results from the LHCb collaboration on semi-leptonic H-meson decays m might provide 
the first indirect hint of new physics beyond the Standard Model (SM). Starting with the 
3cr anomaly in H —)• angular observables [I], several other observables involving 

b —7- sii transitions have produced significant deviations from the SM predictions. The most 
notable is the ratio of —)• to —)• K^e^e~ branching ratios measured as 

Rk = 0-745lQ Qg 2 [2] (deviating from the SM prediction by 2.6a), as in this case the SM 
prediction Rk ~ 1 can be calculated with a very good accuracy. The measured branching 
fraction of Bg —)• [3] is also low compared to the SM prediction. 

These anomalies could well be the result of statistical fluctuations, experimental problems, 
underestimated hadronic uncertainties, or a combination of all three. Nevertheless, it is 
intriguing that many of the discrepancies can be simultaneously explained by assuming new 
physics contributions to the Wilson coefficients of 4-fermion operators with a b- and s-quark 
and 2 leptons [IHIS]. Such contributions can be easily generated in explicit models through 
the exchange of a new Z' gauge boson [T6l - f26] or leptoquarks [IIlETHall. In order to address 
the Rk anomaly, lepton flavor universality has to be broken, which generically also implies 
lepton flavor violation (LFV) [HUl [55H55] . It is tempting to connect these patterns of flavor 
violation to the SM flavor puzzle, i.e. the experimentally observed hierarchical structure of 
Yukawa couplings. Some papers in the literature have attempted to obtain this connection in 
the context of partial compositeness da ESI EH] or leptoquark models jSa- Other works [13 
[22] considered gauged abelian symmetries that are able to reproduce some features of the 
CKM matrix (but not the SM quark mass hierarchies). However, up to now, no model has 
been proposed that directly connects the anomalies to the generation of fermion masses and 
mixings. 

The purpose of this work is to provide a predictive model of this kind. We address the 
anomalies in 6 —?■ sii transitions in the context of a light Z' vector boson, whose couplings 
to fermions are governed by an underlying U{2)f symmetry that explains fermion masses 
and mixings. The original U (2) models proposed in the context of supersymmetry |39| 00] 
have been disfavored by precision measurements in the H-factories m, as they predicted 
the relation Vuhl^ch = which was not borne out experimentally. However, it is 

not difficult to modify this prediction with a more general U{1)f charge assignment, as 
demonstrated in Ref. |42j . 

The model that we are presenting here is essentially a non-supersymmetric version of 
the one in Ref. [l2], in which the dominant source of deviations from the SM is due to the 
tree-level exchange of the Z' gauge boson associated to the U{1)f flavor group. Similarly 
to the supersymmetric model, the couplings of the Z' are approximately U{2)f symmetric, 
and flavor violating effects in the quark sector are suppressed by the small CKM mixing 
angles involving the 3rd generation. In contrast to Ref. |l2], we do not demand that the 
U{1)f charges are compatible with SU{5) grand unification. This generalization gives us 
more freedom in the charged lepton sector to address the observed anomalies in 6 —)• sii 
transitions. Once this is achieved, the parametric freedom in the model is to a large extent 
fixed by matching to the observed quark and charged lepton masses and quark mixing angles. 

In our scenario, the deviations in 6 —sii observables arise from a simple pattern of Z' 
contributions to the 4 relevant Wilson coefficients ■ Namely, the new contributions are 

aligned with the SM one (i.e. approximately left-handed) and controlled by a single parameter 
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(the ratio of the U{1)f gauge coupling and the Z' mass multiplied by the bL-dL mixing angle) 
that sets its overall magnitude. Moreover, they interfere constructively in the electron channel 
and destructively in the muon channel. As a consequence, we predict a simple pattern for the 
relevant b —)• sU amplitudes: in the electron channels the SM predictions are rescaled by a 
factor Te > 1, whereas in the muon channels they are rescaled by a correlated factor < 1. 
The current experimental data on 6 —)■ s£i transitions determine the overall normalization of 
the Z' contribution. This in turn fixes the predictions for other flavor-violating observables up 
to 0(1) coefficients that span the parameter space of our model. Comparing that with existing 
constraints from AF = 2 and LFV observables, we obtain bounds on these 0(1) coefficients. 
The strongest ones come from Bg and kaon mixing, electroweak precision measurement in 
LEP-2, and, especially, from /r-e conversion in nuclei. These bounds disfavor large regions 
of the parameter space, but they nevertheless leave enough room to address the S-meson 
anomalies. The corollary is that our scenario will be decisively tested not only by upcoming 
new data from LHCb, but also from near future tests of LFV in /x —)■ 3e decays and ^-e 
conversion in nuclei. Last but not least, if the Z' boson couples to fermions with electroweak 
strength, it is within the kinematical reach of LHC. 

This paper is organized as follows. In Section 2 we define the setup of the U(2)f flavor 
model and use its predictions for fermion masses and mixings to determine the couplings of 
the Z' gauge boson to fermions. In Section 3 we demonstrate that the resulting contributions 
from tree-level Z' exchange to Wilson coefficients controlling b —)• sU transitions allow one to 
address the H-meson anomalies. In Section 4 we study other constraints on the parameter 
space, and show that electroweak precision tests in LEP-2 and /x-e conversion in nuclei provide 
important constraints. We conclude in Section 5. In Appendix A we provide analytical results 
for the eigenvalues and mixing angles of the quark Yukawa matrices. 

2 The Model 

In this section we define our model with a U(2)f flavor symmetry. We first study its predic¬ 
tions concerning the fermion masses and mixings and demonstrate that the observed patterns 
in the quark and lepton sector can be reproduced. Then we discuss the physics of the Z' 
boson associated to the U(1)f factor of the flavor group. This degree of freedom will be the 
origin of lepton flavor violation in the H-meson sector that we discuss in the next section. 

2.1 Flavor Symmetries 

We first consider an extension of the SM with the global symmetry U(2)f = SU(2)f xU{1)f 
acting in the fermion’s generation space. Here we restrict to the effective description involving 
only SM fields and spurions parametrizing the breaking of SU(2)f x U(1)f- We assume that 
the additional degrees of freedom needed to UV-complete this theory are heavy enough not 
to play role in the low-energy dynamics, i.e. the cutoff-scale A of the effective theory is in the 
multi-TeV range. The first two generations transform as a doublet under SU(2)f, and the 
third generation is an SU(2)f singlet. The U(1)f charges of all fermions are treated as free 
parameters for a while; they will be fixed later to reproduce the observed mass and mixing 
hierarchies. The Higgs field is a total flavor singlet. The breaking of the flavor symmetry is 
described by two scalar spurions: (j) transforming as and y transforming as l_i. These 
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fields acquire the following vacuum expectation values (VEVs): 

W = W = (2.1) 

where we assume <C 1. We also define (j) = which transforms as 2_x^- In Table 

we list the held content and their general transformation properties under the flavor group. In 
the next sections we will specify the U{1)f charges Xf needed to reproduce fermion masses 
and mixings. 



Qa 

Ua 

Da 

La 

Ea 

Qs 

U3 

Ds 

Ts 

E3 

H 

4>a X 

SU{2)f 

2 

2 

2 

2 

2 

1 

1 

1 

1 

1 

1 

2 1 

U{1)f 

A® 

Af 

Af 

Af 

Af 

xg 

^3^ 

^3^ 


V 

0 

A,^ -1 


Table 1: The field content and U{2)f quantum numbers. 


As the fermions are in general charged under U{2)f, Yukawa couplings require additional 
spurion insertions in order to be C/(2)i7-invariant. This leads to non-renormalizable inter¬ 
action suppressed by appropriate powers of A. After inserting the spurion VEVs the cutoff 
dependence drops out, and Yukawa hierarchies arise from powers of the small parameters 
The resulting Yukawa matrices are of the form 


Vf 




,/ ixf+xfi 

,/ \xf+xf-x^\ 


U ,2jXF+Xf+2X^\ 

'^22^<j)^X 

,/ \Xf+Xf+X^\ 


,/ \xj’+xi+x^\ 


hi 


jxf+xi\ 


33'^X 


( 2 . 2 ) 


where P = Q,L, and S = U,D, E. In each entry we omitted terms suppressed by more 
powers of coming from higher-dimensional terms in the effective theory. The absolute 
value appears because only positive powers of y or x* are allowed in the effective Lagrangian. 
Note that, in contrast to the supersymmetric U{2) model in Ref. jl2], there are no holomorphy 
constraints, which leads to a more general Yukawa pattern. 

We move to discussing the consequences of the Yukawa pattern in Eq. (2.2) for the fermion 
masses and mixing. 


2.2 Quark Masses and Mixings 

In the quark sector, we fix = X^ = 0, so that the top Yukawa coupling is not suppressed. 
Furthermore, we impose the following constraints on the charges: 

A^<0, A® + A0>O, Af + A^>0, Vg^YO, Vf + > 0. (2.3) 

With these constraints, we find the following up- and down-quark Yukawa matrices: 



y. u ^11 

0 ^ 


/ 

0 


ud U €32 
^12^12 

0 

\ 

“^12^12 

uu 

^^22^23^32 

^23^23 

, Vd ~ 


hd u 
'''12^12 

^32 

ud u d 
'''22^23^32 

ud ^d 

^^23^23^33 

V 0 

uu 

'^32^32 

^33 / 


V 

0 

32 

u d ^d 
''-32^32 

/j ^ 

'^33^33 

y 
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where we have defined 


vQi -yU 

^12 = 


^23 — ^<t>^X 


^32 — ^<t>^X 


,d = 
^32 — 


^33 = ^X 


(2.5) 


The yii, 2 / 13 , and ysi entries are not exactly zero but one can show they yield subleading 
corrections to quark masses and mixings relatively suppressed at least by e^. Thus, effectively, 
three texture zeros appear in the Yukawa matrix, much as in the supersymmetric models [12] . 
It was pointed out long ago [3^ that the presence of these three texture zeros leads to relations 
among quark masses and mixings that work remarkably well from the phenomenological point 
of view. 


Yukawa matrices of the form in Eq. (2.4) can be diagonalized in a fully analytic way. 


However, it is more convenient to first illustrate the most important points in perturbative 
analysis. Indeed, since mixing angles in the left-handed (LH) quark sector are known to be 
small, one can use them as the small parameter in which the eigenvalues and the remain¬ 
ing mixing angles are expanded. Ignoring 0(1) coefficients, this gives the following rough 
estimates for the eigenvalues and the CKM matrix: 


yt ~ 1, 
Vb ~ £^3 > 
~ e^3 , 


Vc ~ £23^32 ) 


A 


Vu 


'^12'^12 




Y2 


■^23 "^32 


■^23 "^32 
d 

Vd ~ 

^32 

Kb ~ , 

^32 


( 2 . 6 ) 


We thus have 5 small parameters that set the order of magnitude of 8 observable quanti¬ 
ties. Expressed in powers of the Cabibbo angle A ss 0 . 2 , the magnitudes of the parameters 
consistent with experiment is 


£23 ~ £32 


£32 


£33 




-12 




(2.7) 


where the ysjVc hierarchy must be explained by order 1 factors. The parametric size of 
rotation angles and matrices (see Appendix for our conventions) is then given by 


Lu ^ 
12 

^ ^Ru \2 

Lu \ 4 

Sl3 ~ A , 

Ru \ 4 

■S13 ~ A , 

Lu \2 

■*23 ~ A , 

Ru 

*23 

M ^ 
'12 

^ ^Rd \ 

- S12 ~ A, 

Ld \ 3 

■513 ~ A , 

^Rd \ 

Sl3 ~ A , 

Ld \ 2 

^23 ~ A , 

„Rd 

*23 


1 


( 2 . 8 ) 


and 





yf. 


1 A A^^ 
A 1 A2 
,A3 A2 1 



(2.9) 


The constraints in Eq. (2.7) are recovered (up to a mismatch in 6^2 that again is ascribed to 
order one factors) when the U{1)f charges are fixed as 


= 0 , 


= xY = xf = Af = -X^ = 1 , 


( 2 . 10 ) 


and the spurion VEVs are of the order £x ^ £</> ~ A^. One robust conclusion is that, given the 
observed masses and mixings, the C7(1 )f charges in the RH down sector should be universal. 
This has important consequences for phenomenology, as we will discuss later on. 
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We now improve on the above rough estimates, taking into account the 0(1) coefficients. 
For our purpose, it is convenient to express the observables in terms of physical quark masses 
and the unitary rotations that connect the flavor and mass basis. To this end we pick 4 
rotation angles: S 23 , As we show in Appendix, the remaining rotation angles, 

up to percent corrections, can be expressed in terms of these 4 angles and the quark mass 
ratios: 


Lu 

*12 


Ld 

*12 


Ru 


„Rd 



Ijifi 

ms 


Lu 

'’13 


Ld 


Lu Lu 
■'’23 '’12 


Rd 

’23 


_ Ld Ld I _ ^23 

'S23'Sl2 ( J- Rd Ld 

'"23 '’23 ' 'o 


^Ru ^ Ru Lu 
'’13 ~ '’23 '’12 ) 


Rd 

Rd ^ 23 Ld ('9 1 11 

<"23 


Using Eq. (2.11), the CKM elements up to phase factors can be expressed as: 


4 “' 4 “) y m{ cii' 


In the original U (2) models |39( 00] S 23 was taken to be small, ^23 ~ ^cb- From Eq. ( |2.12[ ), 
this leads to the prediction \Vub/Vcb\ ~ \/'mu/mc which deviates from experimental data by 
more than 3a. However, with a large RH 2-3 rotation angle, I/V 2 , the CKM 

angles can be well fit (see also Refs. mills]). The other 3 rotation angles parametrizing the 
model can be small, ~ ~ ~ lUcfel- 


One can explicitly verify that the charge assignment in Eq. (2.10) allows one to fit the 
masses and mixings in quark sector with coefficients that are indeed 0 ( 1 ), and in turn 
check the validity of the above parametrization. For our fit, we take the masses and mixings 
calculated in the SM at the scale 10 TeV 1431. With the Yukawa matrices 


Vu ~ -4.9 • 

V 0 


4.9 -el 0 \ 

3.7-e^ 0.89-6^ 
1.3 • 0.79 / 


and the spurious VEVs 


Vd 


/ 0 3.6 • 0 \ 

-3.6 • e2 -0.62 • 4.9 • 

( 0 -0.20.4 2.5.47 


Ri 0.0040, 


pa 0.035. 


(2.13) 


(2.14) 


one can reproduce the observed masses and mixings within the experimental errors. One can 
also check that the above Yukawa matrices give values for the 2-3 mixing angles 

|4g«1.0'|Kb|, |4“|«1.5'|Kfe|, 14^1 « 2.1-iKbl, |4‘'|«0.60, (2.15) 


Moreover the other rotation angles are in very good agreement with the approximate expres¬ 
sion in Eq. (2.11). We will use Eq. (2.15) and Eq. (2.14) as a reference point for the natural 
values of the mixing angles in phenomenological analyses below. 

In summary, the quark masses and mixing angles can be successfully fit in our U{2)p 
flavor model with 0(1) Yukawa coefficients in the Lagrangian. The remaining freedom can be 
parametrized by four rotation angles. Their precise values depend on the 0(1) coefficients, 
however, barring large cancellations, their order of magnitude is fixed: sg^ ~ 

\Vcb\, and sg* ~ 1 . 


~ sg^ ~ 
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2.3 Charged Lepton Masses 

In the lepton sector we focus on the charged lepton masses, and we ignore here the neutrino 
masse^ Therefore the rotation angles in the charged lepton sector are not constrained by 
phenomenology, which leaves more freedom in the choice of the model parameters. The 
simplest possibility is to take the U{1)f lepton charges to be compatible with SU(5) grand 
unification, that is to say, the same as for the down-type quarks [15]. However, one can show 
that such a choice does not allow to address the lepton non-universality in h —>■ sU transitions, 
which is the primary goal in this paper. Therefore we make a different choice of the U{1)f 
charges: 

Xf = -3-Xf, Xf = 2-Xf, Xf = 4-FXf, (2.16) 

where Xf does not enter into the lepton mass matrix and is left unspecified for the moment. 
This choice leads to the following lepton Yukawa matrix; 


^124 

Ve ~ I 0 1 , 


(2.17) 


where the 2-2 and 3-3 diagonal elements are suppressed by 


44 


and thus can be 


neglected. The consequence is that the muon and tau Yukawa couplings are set by the 23 
and 32 off-diagonal elements. Indeed, diagonalizing Eq. (2.17) yields the Yukawa couplings 

Ufi ~ Ut ~ 


2 /e 


7,e 2 


Using in Eq. (2.14), the correct lepton masses are recovered by fixing three 0(1) coeffi¬ 
cients as hh ~ 0.57, / 1 I 3 ~ 0.29, h ^2 ~ 4.3. The rotation angles are then determined by the 
remaining 0 ( 1 ) coefficients: 



^ 'n3 2 ^ 
" /,e ^X" 
"'23 

5.5 X 10"®/i?3 , 

„Re 

S12 ~ 

^ 4.2 X , 

"32 


^23 " 

-1, 


S23 

0, 



^32 ^<t> 

1.1 X 10"^/l^2> 

Si3 ~ 

_^^_X _ _g 2 X 10“®/l^2 • 

^23 

(2.19) 


In summary, the charged lepton masses can be well reproduced with the U{1)f charge as¬ 
signment in Eq. (2.16). The resulting structure of the Yukawa matrix in Eq. (2.17) leads to a 
large LH mixing between the 2nd and 3rd generation, S 23 ~ 1, and a small RH 2-3 rotation, 
^23 ~ 0 - The remaining freedom is the charge and the three 0 ( 1 ) coefficients hf 2 , ^ 13 ; 
/133 that set the magnitude of the 1-2 and 1-3 mixing angles. We will use this freedom later 
when addressing the 5 —)• s anomalies in a way that avoids phenomenological constraints. 


2.4 Z'-Boson 

We now extend the model by promoting U{1)f to a local symmetry (as in Ref. [42|1. We 
assume that the associated gauge boson is relatively light, with a mass in the TeV range. 

^This is mainly due to simplicity; there are no obvious obstacles to reproduce neutrino masses and mixings 
with Dirac neutrinos upon adding RH neutrinos. 
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Note that the U{1)f symmetry without additional fermions is necessarily anomalous if the 
model explains fermion mass hierarchies. This is due to the relation 


det yuVd ~ ex 


4X5^ +2XY+2XY+2XY +XY+XY 


— 


( 2 . 20 ) 


where C 3 is the anomaly coefficient of the mixed SU{3)‘^U{1)f anomaly. As U{ 1 )f is spon¬ 
taneously broken by the VEVs of cj) and y at a scale v' ~ eA, we assume that the anomaly is 
cancelled by unspecified dynamics (involving new chiral fermions) at the scale A^/y < dvruh 
Since the new gauge boson has a mass given by Mz' = g'v', it can easily be the lightest 
new degree of freedom when g' is sufficiently small. We therefore ignore the additional heavy 
dynamics and concentrate on the effects of the Z' gauge boson. 

In the flavor basis, Z' couples to each fermion proportionally to its U{1)f charge Xf, 

C D g'Z'^ + X^u} + X^Dja^^Di + X^Llaf^U + Xf . (2.21) 

We have fixed these charges (except for Xf) to fit the observed fermion mass hierarchies. 
That fit also determines the unitary rotations that connect the flavor and the mass basis. 
Therefore, flavor non-universal effects mediated by Z' are predicted in our model, up to an 
overall normalization determined by the Z' mass and gauge coupling, and up to the freedom 
of choosing Xf and order one Yukawa factors. In particular, the SU{2)f structure for the 
first two generations implies that flavor changing effects are entirely determined by the 3rd 
row of the rotation matrices. In the mass basis, the Z' couplings take the form 




( 2 . 22 ) 




= Xf 




Sij 


Xf 

Xf - Xj 

xf 




V = Q,L, 

S = U,D,F. 


(2.23) 

(2.24) 


Note that flavor-violating couplings are proportional to the charge difference Xi — X 3 . As a 


consequence, with the charge assignments in Eq. (2.10) there is no flavor violation in the RH 
down sector: = 6 ij . Eor the LH down quarks we hnd 


A didj 

L 


1 

nLd r^Ld 

„Ld \ 
*13 ) 

, / 1 

A3 

A3 

^L/d ^Ld 
*13 *23 

{cki? 

ski 


1 

A2 

„Ld 

*13 

„Ld 

*23 

(ski?) 

' V A3 

A2 

A4 


(2.25) 


Since sf '3 <C S 23 * (see Eq. (2.11)), the largest flavor violating effect of Z' is in 6 —)■ s quark 
transitions. This will be handy for addressing the recent R-meson anomalies, as we will 
discuss in the next section. Eor LH and RH up quarks we find 


A UiUj 

^x 


1 

r>XU aXU 

^13 *^23 

„Xu \ 
*13 i 

, / 1 

A® 

A^ 

„Xu „Xu 
*13 *23 

(4“)^ 

*i“ 

~ A« 

1 

A2 

„Xu 

*13 

„Xu 

*23 


' V A^ 

A2 

A^ 


(2.26) 


where X = {L, R}. 
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In the charged lepton sector the Z' couplings depend on the charge Xf and to good 
approximation only on the 1-3 rotation angles. 


A eiEj 


-3 - Xf 

r.„Le 

0613 


0 \ 


/ xf 

0 

46f| 

2 - Xf 0 


0 

xf 

0 

0 -3-Xf / 

\ 

0 

4 + Xf 


(2.27) 

The diagonal muon Z' coupling is different from the electron and tau one at leading order, 
which is due to the 2-3 inversion in the LH sector. This feature of our model will allow us later 
to address the anomalies in i?-meson decays involving muons and electrons. In the RH sector 
electrons and muons have approximately the same coupling to Z' and the dominant flavor 
non-universal effects must involve the tau lepton. The largest flavor violating effects occur in 
the LH fi-e and RH r-e transitions. Note that the rotation angles setting the magnitude of 
these lepton-flavor-violating effects are fixed up to 0(1) factors and expected to be tiny, see 


Eq. (2.19). As a result, lepton flavor violation in our model is suppressed at least by a factor 
of order e^/e<^ ps 5 x 10“^ as compared to violation of lepton flavor universality. 

3 Phenomenology of 6 —> s£i Transitions 

We now turn to the predictions for b —)• sii transitions. Our main goal is to address the 
recently observed violation of lepton flavor universality in R-meson decays [2] . In our model, 
this anomaly is due to the exchange of a U(1)f Z' boson with mass in the multi-TeV range. 

Low-energy observables are controlled by the 4-fermion effective operators that arise from 
integrating out the Z' at tree level, 
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/2 






. fifj 


,fkfl f, = 


, fkfl 


(3.1) 


The 4-fermion operators in Eq. (3.1) include the ones relevant for B —>• Kil decays which are 


customarily parametrized by the following effective Hamiltonian (see e.g. Ref. 


Heff 3 [s^^PLb] [cfiri' + 


+ h.c. 


(3.2) 


Note that the analogous 4-fermion operators with RH quarks are not generated in our model 
(Cf = Cfo' = 0), as a consequence of the universal U(1)f charge assignment in the RH down 


sector. Matching Eq. (3.1) and Eq. (3.2), the Wilson coefficients are given by 

J2 , , 


r'U _ 
Ug — - 


Fg- 


V2Ml,GFVtbV^*a 


Af ( Af' + Ai 


riW _ 
'-^10 — “ 


Fg 


V2Ml,GFVtiyt*sa 


Af f A^/ - A"' 


L 

(3.3) 


Eocusing for now on the lepton flavor conserving operators with electrons or muons, we have 


G ee 

9 


0.19A:. 


Cr ~ -0.13A:, 


Gll « -0.19 (1 + 2/3Af) k , 
Ri 0.13 (1 - Xf) k , 


(3.4) 
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where we defined the parameter k as 


/20TeV 

\Mz'/g' 



(3.5) 


and we used the numerical values |Vts| \Vcb\ ~ 0.041, a{mi)) « 1/133. 

The explicit expressions for the Wilson coefficients in Eq. ( |3.4[ ) imply that large correc¬ 
tions to i?-meson decays involving muons are correlated with comparable corrections to the 
analogous observables involving electrons. While new physics in the muonic sector alone gives 
the most economical explanation of the LHCb anomalies (including Rx |12l I14| i. it has been 
emphasized that large corrections in electron channels are not only allowed but could also 
(slightly) improve the goodness of the global fit jl3| . 

We can now identify the parameter space of our model where the measurements of semilep- 
tonic b —)• sii transitions with i = e, g are best reproduced. As can be seen from Eq. (3.4), 
these observables depend on just 2 parameters: k defined in Eq. (3.5), and the U[1)f charge 
Xf. For the moment, we treat them as free parameters, various precision constraints will 
be discussed in the next section. In order to find the best ht region for k and Xf, we use 
the result of Ref. [13]. The authors provide the results of the fit in the 2D planes Cg^), 
{Cfg, C^q), {Cg^, C^g) und (C|®, Cfg). Ignoong possible correlations in the full 4D likelihood, 
we identify the allowed range for k for discrete values of Xf by requiring to simultaneously 
remain inside the 68% or 95% confidence level regions in every 2D plane. We obtain: 


xf 

Icr Region 

2 c7 Region 

-3 

- 

k G [0.5, 2.1] 

-2 

- 

k G [0.5, 3.1] 

-1 

- 

k G [0.6, 5.1] 

0 

k G [2.7,4.2] 

k G [0.8, 6.6] 

1 

- 

k G [1.2,4.9] 

2 

- 
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This simplified analysis suggests that, while a reasonable ht to the b —)• sii data is possible 
for a range of Xf, the best case scenario is X^ = 0. In the rest of this paper we focus 
on that particular choice. In this case, the new physics contributions mediated by Z' are 
purely left-handed, Cg® = —Cfg, Cg^ = —Clfg, and we can derive constraints on k using the 
2 -parameter ht of Ref. [T3| for precisely this case. This way, we hnd that the Icr conhdence 
interval is k € [1.9,4.9] and the 2 ct one is k € [0.2, 6.5]. The allowed region of the parameters 
Cg^-Cg^ parameter space overlaid with the prediction of our model is displayed in Fig. 

We turn to discussing predictions of our model. The case with Xf = 0 is particularly 
simple because the SM contributions to the effective Hamiltonian in Eq. (3.2) are also purely 
left-handed, Cg^ ~ with Cg^{mb) ~ 4.2. Therefore the new physics contributions 

interfere constructively with the SM in the e-channel and destructively in the /r-channel, 
resulting in a simple rescaling of R-meson decay rates by the factors r^ and that are the 
same for all b —see and b —)• s///i processed 


re « (10.044A:)^ r,, ^ {1 - 0.029k f . 


(3.6) 


2 


This a good approximation in the limit where contributions from fully-hadronic operators can be neglected. 
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Figure 1: Global fit results from [13], the blue (light blue) domain corresponds to the Icr {2a) 
region. The red line is our model prediction for Xi = 0 and varying k. 


For k G [1.9,4.9] we thus predict an enhancement in all the electron channels by G 
[1.17,1.48], and a suppression in all the muon channels by G [0.73,0.89]. As a refer¬ 
ence, in Table we show the measured values and the SM predictions for various b —?■ sii 
observables. It is remarkable that most observables in the muon channel shows a deficit com¬ 
pared to the SM predictions, while the ones in the electron channels show some (albeit not 
statistically significant) enhancement. 

Lepton flavor universality is often tested by measuring ratios of branching fractions of 
semileptonic i?-meson decays. In our model we have 


B{B^X^'^fj, ) /I —0.029A:\^ 

B{B Xe+e-) ^ ^ \1 + O.OUk) ’ 


where X = K,K*,(j),Xs. The interval k G [1.9,4.9] corresponds Rx G [0.50,0.76], which 
should be compared to the LHCb measurement Rk = 0.745^q Qg 2 . Future improvements in 
the precision of Rk and other measnrements will be crucial for testing our model, since we 
predict a rather low value for these observable. This is actually supported by measurements 
of inclusive B —> XgM decay ratios from BaBar {Rxs = 0-58 ± 0.19) |17| and Belle {Rxs = 
0.42±0.25) |51j . although with large errors. Another test of lepton-non-universality is provided 
by double ratios |52| such as Rk*/Rk- In our model, as in any scenario with Cg^io = 0, all 
these double ratios are predicted to be equal to one. 

Finally, we comment on the predictions concerning LFV B-meson decays. The Wilson 
coefficients of 4-fermion operators mediating these decays are suppressed by additional powers 
of small parameters, e.g. 

<^r ~ -Cfo ~ -0-3 sfl k , (3.8) 
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Observable 

[GeV^j 

SM prediction 

Measurement 

Ratio 

10«(GeV)^ X B{B+ X^e+e") 
10®(GeV)2 X B{B+ X^e+e") 

[1.0,6.0] 
[14.2,25.0] 

1.73 ±0.12 

0.20 ±0.06 

1.93 ±0.55 gT] 
0.56 ±0.19 gT] 

1.12 ±0.33 

2.80 ±1.27 

10^ X B{Bs —^ 

109(GeV)2(^)(H0 ^ XV+M“) 
109(GeV)2(^)(H0 ^ XV+M") 
109(GeV)2(|4)(H+ ^ K+^i+ijl-) 
109(GeV)2(^)(H+ ^ 
109(GeV)2(^)(H+ ^ 
109(GeV)2(^)(H+ ^ X>+^-) 
106 (GeV )2 iB{B+ X,/x+/i-) 

106(GeV)2 X B{B+ X,/x+/i-) 

[1.1,6.0] 
[15.0,22.0] 
[1.1,6.0] 
[15.0,22.0] 
[1.1,6.0] 
[15.0,19.0] 
[1.0,6.0] 
[14.2,25.0] 

3.54 ±0.27 

31.7 ±9.4 

13.6 ±2.0 

34.8 ± 10.3 

14.8 ±2.0 

50.5 ±28.6 

61.5 ±34.8 

1.66 ±0.12 

0.24 ±0.07 

2.9±0.7 gSlIlHI 
18.7 ±3.6 [50| 

9.5 ± 1.7 |50| 
24.2 ± 1.4 [50] 
12.1 ±0.7 gUj 
36.6 ±8.7 go] 

39.5 ±8.5 go] 
0.66 ±0.88 g7] 
0.60 ±0.31 g7] 

0.8 ±0.2 

0.59 ±0.21 

0.70 ±0.16 

0.70 ±0.21 

0.82 ±0.12 

0.72 ±0.45 

0.64 ±0.39 

0.40 ±0.53 
2.50 ±1.48 


Table 2: The measured values and the SM predictions for various b —)■ s£i observables used 
in the fit of Ref. [13] . 


Decay 

Branching ratio 

Ref. 

Type 


(4.29 ±0.22) X 10"^ 

m 

LEG 

R+ ^ X+/i±r=F 

< 4.8 X 10“^ 

g3| 

LEV 

B+ K+e^T^ 

< 3.0 X 10-5 

m 

LEV 


< 9.1 X 10"® 

m 

LEV 

Bs 

(2.9 ±0.7) X 10-^^ 

mm 

LEG 

Bs e^ijF 

< 1.1 X 10-® 

m 

LEV 


Table 3: Experimental constraints on the branching fraction of several lepton flavor conserv¬ 
ing (LFC) and lepton flavor violating (LEV) R-meson decays. 


From Eq. (2.19), typical values of the lepton mixing angle are ®13 ~ 10 ^-10 This implies 
the rate of lepton flavor violating decays is suppressed by at least 6 orders of magnitude 
compared to the lepton flavor conserving ones. Given the present sensitivity summarized in 
Table this will not be observable in the near future. LEV decays involving tau leptons are 
even more suppressed. 


4 Constraints from Flavor Violation and Direct Searches 

In this section we discuss constraints on the parameters of our model from AF = 2 flavor 
transitions, LEV decays of leptons, LEP-1 and LEP-2 electroweak precision observables, and 
direct Z' searches at the LHC. We will show that it is possible to address the observed 
violation of lepton-flavor universality in H-meson decays without violating these constraints. 
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4.1 Constraints from AF = 2 Observables 


The Z' boson exchange generates 4-quark operators mediating AF = 2 transitions. In the 
notation of e.g. Ref. [55] these are denoted Oi (with 4 LH quarks), Oi (with 4 RH quarks), 
and O 5 (with 2 LH and 2 RH quarks). Since = 0, for the down-type quark only Oi is 

generated. Their (in general complex) Wilson coefficients are given by 


Cf = 


g'2^sd^sd 


2M|, 


Cf'* = 


g'^A¥A¥ 


2M|, 

and numerically one has (using md/mg ~ 0.05) 

1.8 X 10"^° 


Cf ^ = 


2M% 


rB, 1-1 X 10 

^ TeV^ 


Cf = 
Vl^cfeiy 


-7 


TeV" 


C23 R ^ *-"1 — 


„Ld 

■^23 

l^cbl 


^2i k , 

2.1 X 10 

TeV^ 


Bs 


-6 


Vl^cfeiy 


(4.1) 


(4.2) 


(4.3) 


where k is defined in Eq. (3.5) and it needs to be 0(1) for the model to address the B- 
meson anomalies. We have also approximated sf'g ~ ”'* 23 'Si '2 > which slightly overestimates 
the Wilson coefficients, see Eq. (2.11). These expressions have to be compared to the bounds 
from iL-mixing taken from Ref. [55| (Im) and Ref. [56] (Re), and the bounds from R-mixing 
taken from Ref. m- 


ImOf < 

l<^fl < 


3.4 X 10-9 
TeV^ ’ 

1.4 X 10-® 
TeV^ ’ 


Re of < 

I of “I < 


9.6 X 10 
TeV^ 
1.8 X 10 
TeV^ 



(4.4) 

(4.5) 


This shows that for S 23 ~ Reft and k in the experimentally preferred range k G [1.9,4.9] the 
bounds from K, B and Bg mixing are satisfied, even for an 0(1) phase in ImOf. 

Turning to the up sector, 4-fermion operators with both LH and RH fermions are gener¬ 
ated, as a result of a non-universal U{1)f charge assignment. In particular, for the AO = 2 
operators we have (with mu/rric ~ 0 . 002 ) 


-iD 


■-iD 


-,D 


7.1 X 10-^2 


TeV^ ' 

URcfel 

7.1 X 10-12 

/ „Ru 
f '’23 

TeV2 ' 


2.8 X 10-11 

( aii 

TeV2 ' 

URcfel 


^23 


'’23 

„Ru 

^23 


m) \s^iJ 


(4.6) 


Hence, D-meson mixing is further suppressed compared to K-meson mixing by the small mass 
ratio Given the bounds from Ref. |58j (Im) and Ref. [56] (Abs): 


ImOf < 


0.9 X 10-® 
TeV^ ’ 


ImOf < 
ImOg^ < 


0.9 X 10“® 
TeV^ ’ 
0.4 X 10-® 
TeV^ ’ 


[Of 

[Of 


7.2 X 10-'^ 

TeV^ ’ 
7.2 X 10“^ 
TeV^ ’ 
4.8 X 10-^ 
TeV^ ’ 


(4.7) 
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Decay {ij)ikl) A^'=^'|/(V2 


Ks —)• e+e 
Kl —)■ e'''e“ 

Ks 

Kl -> 

—>■ 7r+e'’'e“ 
—>■ 7r°e'''e“ 
—)• 7r+/i+//“ 

S'*" —>■ 7r'’'e'’“e“ 
B~^ —)• 7r+A/''’A^~ 


(21)(11) 

< 1.0 

(21)(11) 

< 2.7 X 10 

(21)(22) 

< 5.1 X 10 

(21)(22) 

< 3.6 X 10 

(21)(11) 

< 6.7 X 10 

(21)(11) 

< 1.6 X 10 

(21)(22) 

< 5.3 X 10 

(31)(22) 

< 3.9 X 10 

(31)(11) 

< 2.8 X 10 

(31)(22) 

< 2.3 X 10 


Table 4: Upper bounds on Wilson coefficients from leptonic and semi-leptonic K and B 
decays with s —> d and b ^ d transitions. 


there are no further bounds on our model from D—mixing. 


4.2 Semileptonic decays in 6 —A d and s —)■ d transitions 

We now turn to semileptonic decays involving 6 —)• d and s —)• d transitions with electrons 
or muons in the final state. The relevant Wilson coefficients for the associated 4-fermion 
operators for 6 —?■ d transitions are given by 


r(M)(ee) ^ ^ 3.4 X 10-5 

2M|, TeV2 ’ 

C(m)M ^ ^ 2.3 X 10-5 ^ 

2 M|, TeV^ 

and for the s —>■ d transitions by 

risd)iee) ^ g'^KfKf ^ 1.4 X 10-6 

2M|, TeV2 WcblJ 

risd)M = ^ 9.4 X 10-^ (s^\ 

2M|, TeV^ V|4^cb|y ’ 


(4.8) 


(4.9) 


where we approximated A^*^ —s^i^/md/ms and AJ^ (s 23 )^ -^/md/ms- The upper bounds 
are summarized in the Table adapted from the case of composite leptoquarks [28] . From 
this table, it is easy to verify that for the experimentally preferred range k G [1.9,4.9] all the 
bounds are satisfied. 


4.3 Lepton flavor violation 


From Eq. (2.27), the largest LFV Z' couplings are the ones to LH muons and electrons. These 
are constrained by several precise measurements of LFV p —)• e transitions. First, we have 
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the /i —>■ 3e decay with the branching fraction 


/4 4 
g V 


BR(f, ^ 3e) = (2|A"P + 1^511") 


“■z- 
3.2 X 10 




„Le 

^13 


(\Vcb\y 

V sM ) 


(4.10) 


1.1 X 10-4 

where v = 246 GeV. This should be compared with the experimental limit from Ref. [59] : 

BR(^ ^ 3e) < 1.0 X 10"^2 _ (411) 


This limit can be violated for larger values of the parameter k. In particular, for = 0 
and k in the experimentally favored range k G [1.9,4.9] we get the constraint on the mixing 
angles: 


V 4 v Vi-ixio 


„Le 

^13 


< [3.0, 1.2]. 


(4.12) 


This can be satisfied for ~ \^ch\ if the 0(1) Yukawa coupling controlling is < 1 - 

A stronger constraint on LFV comes from g-e conversion in nuclei. Borrowing the formulas 
e.g from Refs. |60H62j . for Xf = 0 the conversion rate in gold, titanium and aluminium nuclei 
is given by: 


CR(m ^ .. Au) » 4.0 X 10- (!1|1)' ( , 

CK(, ^ .,Ti) = 3.8 X 10-- (!i|l)' 

CR(,. ^ e, Al) = 1.0 X 10-- (!i|l)' (4.13) 

This should be compared with the bounds from Ref. |63[ : 

CR(/x ^ e,Au) < 7.0 X 10■^^ CR(/4 ^ e, Ti) < 4.3 x 10“^^ (4.14) 


For the parameter k in the range favored by the R-meson anomalies, k G [1.9.4.9], this leads 
to the constraint on a combination of mixing angles in our model 


) Vl-1 X 10 


„Le 

^13 


< [0.22,0.08]. 


(4.15) 


Formally, 5^3 is a free parameter, therefore Eq. (4.15) can always be satisfied with an appro¬ 
priate choice of the lepton Yukawa couplings. However, our philosophy is to explain the flavor 


IS sfo® ~ 


’13 


hierarchies with all Yukawa couplings in Eq. (2.2) being 0(1), in which case the natural value 
10-4. In this respect, Eq. (4.15) forces us into a less natural corner of the parameter 


space and suggest a value of k close to the lower la boundary. We note that the experimental 
sensitivity to the /i-e conversion rate is expected to improve by many orders of magnitude in 
the near future |65H67j . In case of a null result, our model will no longer be an attractive 
solution to the R-meson anomalies. 
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4.4 Electroweak precision tests 


Integrating out Z' induces lepton-number conserving 4-fermion operators which can be con¬ 
strained by electroweak precision tests. Here, we focus on the 4-lepton operators which give 
the strongest bounds due to large U{l)p charges of leptons. At leading order, these do 
not affect Z-pole observables measured in LEP-1 and SLC, but they can be constrained by 
off-Z-pole fermion scattering in LEP-2. We parametrize these operators as 




E 

£&e,ii,T 


For Xf = 0 the non-zero Wilson coefficients are 
[cLLjee = “ 


2M2 


(Af = - 6.8 X 






Xd 

’23 


[cLL]ev — ^[cLL]eei 


[cLL]eT = ‘2'[cLL]e 


8 . 


[clrIst = --[cLL]e 


(4.16) 


(4.17) 


(4.18) 


Note that the sign of each contribution is hxed, in particular the contribution to [c^^Jee is 
always negative in our model. We calculated the impact of these operators on the LEP-2 
observables quoted in Ref. [Mj- We used the total cross section and asymmetries of e''"e“ —>■ 
measured at the center-of-mass energies ^/s G [130,207] GeV, as well as the 
differential cross-sections of e'''e“ —)• e"''e“ at ^/s G [189, 207] GeV. This way we obtain the 
95% CL constraint: 


\ ^L/d 1 
V *23 / 


< 1 . 1 . 


(4.19) 


For k in the experimentally favored range k G [1.9,4.9], Eq. (4.19) requires somewhat larger 
values of of order 2-4 |I4fe| • This leads to some tension with the bound from CP violation 
in kaon mixing in Eq. ( |4.4[ ), assuming 0(1) phases entering . Much as the LEV bound, 
these constraints point to rather low k k, 2. 

We also comment on the corrections to lepton flavor conserving muon decays. Loops with 
a Z' boson result in the following 1-loop correction to the /i —?■ ev^Ve decay width 

r(/i —)■ evv) 


r(^ 


env)sM 


= 1 - e, 


(4.20) 


where 


/2 


e Ri — 


^9 A ee 

47r2 


Ml 


= 1.5 X 


W)(S) ■ 


The muon lifetime measurement does not constrain new physics by itself, because it is used to 
extract the SM input parameter Gp (equivalently, the Higgs VEV v). However, indirectly, new 
physics contributions to Gp shift other observables (for example Z-pole asymmetries, 
etc.) away from the SM predictions. To estimate the resulting constraints, we note that the 
effect in Eq. (4.20) is equivalent to introducing the 4-lepton operator 

with the Wilson coefficient [cll\i 221 = £• The constraint on this Wilson coefficient from the 
Z-pole observables can be read off using the global likelihood function quoted in Ref. [69]. 


If only this one operator affects the Z-pole observables, the constraint reads —0.8 x 10“^ < 
[ci,l]i221 < 2 X 10“^ at 95% CL. The resulting constraints on the parameter of our model are 
weaker than the ones from off-Z-pole measurements in LEP-2. 
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Figure 2: For = 0, the region of the Mz'-g' plane of our model excluded by resonance 
searches at the LHC (red). We also show the indirect constraints from the 2-fermion production 
in LEP-2 (black mesh). The green regions correspond to s^i = 211461 and the parameter k in 
Eq. (3.5) in the range favored by the i?-meson anomalies at 1 cr fc S (1.9,4.9) (darker) and at 


2 a k G (0,2,6.5) (lighter). 


4.5 Z' searches in colliders 


Finally, the parameter space of our model is constrained by direct searches for resonances 
in colliders. Since addressing the i?-meson anomalies requires Mz'jg' ~ 20 TeV, the Z' 
boson predicted by our model is within the kinematic reach of LHC for g' of electroweak 
strength or smaller. Note that the direct searches probe separately the Z' mass and coupling 
constant, unlike all previously discussed observables that depended on these parameters only 
via the combination Mz’jg'■ Given the charge assignments in Eq. (2.10) and Eq. (2.16), the 
branching fraction of Z' into dilepton hnal states is significant. In particular, for Xf = 0, we 
have 

Br(Z' -G ee) ^ 14%, Br(Z' ^ /i^) sa 6%, (4.22) 


and the strongest constraints are expected from the di-electron channel. In Eig.j^we plot the 
constraints in the Mz'-g' plane based on the CMS search for di-electron resonances in the 
LHC at ^/s = 8 TeV [TOj. These constraints imply Mz' > 3 TeV and g' > 0.1 in the region 
of the parameter space favored by the H-meson anomalies. Note that the direct limits are 
complementary to the indirect ones from LEP-2. The latter would allow us to address the 
H-meson anomalies with a light (m^/ < 2 TeV) and very weakly coupled Z'] such possibility 
is however excluded by the resonance searches. 
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5 Summary and Conclusions 


In this work we have addressed recent anomalies in 6 —)• sti transitions. This was achieved 
by a tree-level exchange of a light Z' vector boson, whose couplings to the SM fermions 
are governed by an underlying U{2)p flavor symmetry that explains fermion masses and 
mixings. While the U (2 )p quantum numbers of quarks are determined by quark masses and 
CKM angles, there is more freedom in the charged lepton sector. However, requiring that 
the b —7- sU anomalies are fit within lu essentially selects unique U{2)p quantum numbers 
for charged leptons. The only free parameters are then the Z' mass, the associated gauge 
coupling and a handful of 0(1} parameters controlling flavor violation. As a result of the 
U{2)p symmetry structure, the magnitude of the latter is set by rotation angles involving the 
third generation, which implies that flavor violation in the quark sector is strongly suppressed 
by small CKM angles. 

In our model, Z' couples only to LH electrons and muons and has universal couplings 
to RH down quarks. Therefore, the Wilson coefficients relevant for the b —)• sii transitions 
satisfy and are therefore approximately aligned with the SM contribution. 

Furthermore, the ratio Cg^fCg^ is fixed by the U(2)p quantum numbers. As a result, b —)> sU 
transitions are governed by a single parameter k 


k = 


f20TeV 



(5.23) 


controlling the overall magnitude of the Z' contributions. Interference of the new contributions 
with the SM is constructive in the electron channel and destructive in the muon channel. This 
leads to the prediction that the semi-leptonic R-meson decay rates are rescaled by factors rg 
and that are common for all b —>■ see and b —)• sgfj, processes, respectively. We find 

( 1 — 0 029A:\ ^ 

TToMikJ • 


Using the fits to all available R-meson data from the previous literature, we determined the 
Icr confidence interval for k: 


fe G [1.9,4.9] , (5.25) 

and thus we predict the rescaling factors in the electron and muon channels, along with their 
ratio Rx 


re G [1.17,1.48] 


G [0.73, O.f 


Rx G [0.50,0.76]. 


(5.26) 


The parameter k also controls flavor-violating processes, but only in conjunction with quark 
and lepton rotation angles which are fixed up to 0(1) parameters. The strongest constraints 
on these parameters arise from /r-e conversion in nuclei. They can be satisfied if the mixing 
angle sfl is somewhat suppressed by these 0(1) parameters, and/or the mixing angle is 
somewhat enhanced. These considerations also favor smaller /c, close to the lower limit of the 
Icr confidence interval. Electroweak precision test from LEP-2 lead to similar conclusions. 
Explicitly, all bounds are satisfied for k k, 2 and the reference values in Eq. (2.15), for an 
0(1) coefficient < 0.5 in sff. 

All in all, there exist regions of the parameter space of our model where the R-meson 
anomalies are explained and other experimental constraints are satisfied. The final verdict 
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will be provided by further of lepton flavor universality in LHCb and B-factories, as well as 
by future experiments looking for fi — e conversion in nuclei. Apart from indirect searches, 
the Z' boson is likely within the reach of the LHC run-2, and should first show up in the 
di-electron channel, as a result of its large coupling to electrons. 
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A Exact Quark Rotation Angles 


In the flavor basis the Yukawa terms read 

C C +Q^yuUH + QvdDH + L^yeEH. (A.l) 

where each yf is a 3 x 3 matrix with in general complex elements. Here, for simplicity, we 
assume all entries are real, as this will allows us to obtain compact formulas for the eigenvalues 
and the mixing angles. Going to the mass eigenstate basis involves unitary rotations defined 
by 

yu = Vtyi“‘HVlif, yi = Viyf‘>{V‘^f, Ve = Vlyf‘HVf.f, (A,2) 

We parametrize the rotations as 


where (I = L, R) and 


Vlo = -s 


Vl = 


V/3 = 


Vr = EjUif Uj, 


VO 0 1/ V-4 0 4, 

We are interested in the Yukawa matrix of the form 


-2/12 2/22 2/23 
0 2/32 2/33 


yk = 


ydiag = 0 y 2 0 


2/1 < 2/2 < 1 • 


Plugging this form into the eigenvalue equations, from the 1-1, 1-3, 3-1, 1-2, and 2-1 entries 
we get 4 independent equations which are 

-I- c{|sf2-Sl3^32 4- cf'3sf'2'Sf3£23 “I" 'Sl3'Sf3e33 “t" Ci2Cl2Cl3c{|2/l = 0 , (A.6) 

Cl3cf3^12^23 -I- cf3Si3e33 — cf'3sf'2'Sl2'Sf3^22 “ ■sf2S 13'513^32 “ Ci2cf2Ci3Si32/l = 0 , (A.7) 

Cl3cf3^f2^32 — cf3Si2'Sl2'513^22 + Ci3sf3e33 — Si2Sl3'Sl3e23 “ Ci2cf2cf3Si32/l = 0 , (A.8) 

^12^13512^22 + Cl2Ci3Si2f22 + Ci25i3e32 + Ci2'Sl3^23 ~ (^12^13^12yi ~ Cl2Cl3'®122/l = 0 , (A.9) 
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where we have defined 


£22 — S 23 S 2 S + C23C^3y2 , 

£32 = C23S23 — £^ 35232/2 , 


£23 — C^'S 23 “ ^238231/2 , (A. 10 ) 

£33 = C23C§3 + 823823112 , (A. 11 ) 


with 

£22£33 

From the first three equations one finds 


£23£32 — 2/2 • 


.R 

*12 


Cf3£32 


and 


L R I ^12^122/2 

C 12 C 122/1 +- 

£33 


*12 


= 0 . 


*13£33 

^^3^23 


Using this last equation, one gets from the 4 th equation the solution 

£33 [(cf' 3 )^£ 232 /i + (*13)^ (-£i32/i + 2 '^^ 


- 

£13 — 


Ci3e232/i2/2 


(A. 12 ) 


(A. 13 ) 


(A. 14 ) 


(A. 15 ) 


Using this solution and Eqs. (A. 13 ) in Eq. (A. 14 ), one finally can solve for (cl'3)^ in terms 
of Cij and yi^2- Although the exact solutions are not very complicated, one can approximate 
these expressions using yi <C 2/2 1 to get the final solution (with a consistent choice of 

signs): 


^ 

(£232/1 


13-y 

£332/2 

£13 ^ 


/ ei22/i 

„R ~ 

13-y 

£332/2 

£13 - 

R 

*13£33 

„R 

12 — “ 

cf3£32 

£12 — 

L 

*13 £33 

L 

12 — “ 

£13^23 

£12 — 


' 1 - 


£ 232/1 

£332/2 


' 1 - 


ek/i 

£332/2 


cfa = sign(-e23£32)\/l - { 8 ^ 2 )“^ 


(A. 16 ) 

(A. 17 ) 

(A. 18 ) 

(A. 19 ) 
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